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Abstract. Preparing molecules at rest and in a highly pure quantum state is a 
long standing dream in chemistry and physics, so far achieved only for a select set of 
molecules in dedicated experimental setups. Here, a quantum-limited combination of 
mass spectrometry and Raman spectroscopy is proposed that should be applicable to 
a wide range of molecular ions. Excitation of electrons in the molecule followed by 
uncontrolled decay and branching into several lower energy states is avoided. Instead, 
the molecule is always confined to rotational states within the electronic and vibrational 
ground-state manifold, while a co-trapped atomic ion provides efficient entropy removal 
and allows for extraction of information on the molecule. The outlined techniques 
might enable preparation, manipulation and measurement of a large multitude of 
molecular ion species with the same instrument, with applications including, but not 
limited to, precise determination of molecular properties and fundamental tests of 
physics. 



1. Introduction 

In the last 30 years, laser cooling[ll [2] and laser spectroscopy has led to spectacular 
developments in atomic physics. The electromagnetic vacuum of infrared to ultraviolet 
light provides a near perfect entropy sink and allows us to optically pump and cool 
atoms, if sufficiently closed cycles of electronic excitation and subsequent decay are re- 
alized. Entropy can be removed from the electronic state and the translational motion 
to where these degrees of freedom are in a highly pure quantum state. Such quantum 
state preparation is available for a sizeable number of neutral and charged atoms, but 
in most molecules the decay of excited electrons branches into a multitude of levels, and 
this prevents efficient entropy removal in all but a few cases. These exceptions include 
cooling molecules by buffer gas^, exploiting fortuitous level structures that allow for 
closed cycles with a manageable number of laser sources |1H6], selective ionization and 
charge transfer [7j or assembling molecules from ultracold atoms [5HTT]. In all of these 
cases, the molecular states are destroyed, or the molecules leave the interaction region 
during the process of their detection. Methods to prepare ensembles of specific molecules 
in their absolute ground state by pulsed adiabatic passage [T2HT1] or by collisions with 
a laser cooled buffer gasp^ have also been envisioned. 
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2. Translational cooling of molecular ions 

In the method proposed here, the molecular ion of interest is co-trapped with an 
atomic ion that can be precisely manipulated by laser and microwave fields [16]. Due 
to the Coulomb-coupling between the co-trapped ions, laser-cooling of the atomic 
ion will sympathetically cool the translational motion of the molecular ion as well. 
Typically, excited electronic and vibrational states of the molecular ion will decay to 
their ground states within nanoseconds and milliseconds, respectively, and the rotational 
state distribution will equilibrate to the black-body temperature of the trap within a 
few seconds E]- After this initial equilibration, the molecule ideally never leaves the 
electronic and vibrational ground-state manifold. The mass of the molecular ion can 
be inferred quite precisely by determining the normal-mode frequencies of motion of 
the co-trapped ions [171 [18]. This can be accomplished by sweeping the frequency of 
an electronic drive signal that is applied to a trap electrode [T7|. If the drive frequency 
coincides with a normal-mode frequency, the crystal gets excited and the laser-induced 
fiuorescence changes. If the mass of the atomic ion is known and it has frequency Ua 
in the potential well, the two- ion crystal with mass ratio fi = rrim/ma > 1 will exhibit 
two axial modes, with frequencies [T9| 




(1) 



and eigenvectors 



v± = 1^ ^ ^ { 1 - /i T + l} , (2) 

I- fiT a/zU^ - ^ + 1) + 1 



where the first coordinate is of the the lighter mass atomic ion. If the molecule mass 
rrim is smaller than rria similar relations can be derived. From the measured frequencies 
Ua and u±, which can be determined with better than 10~^ relative accuracy p^, we 
can infer m^: 

mm = rria^— — / o 2 " W 

For /i ^ 1, goes to \/2 to a while u;_ goes to zero. In the same limit the normal- 
mode amplitudes of the ions go to {—1,0} in the u+- mode and to {1/a/5, 2/-\/5} for 
the a;_-mode. Cooling of both modes will become more of a challenge with larger mass 
differences: For the a;+-mode, the participation of the heavier ion tends to vanish, while 
the a;_-mode will go to low frequencies, where resolving sidebands becomes harder and 
ambient heating mechanisms are more detrimental [TB]. 

For not too different masses, it should be possible to cool the normal modes of 
motion of both ions to the ground state, analogous to ground state cooling demonstrated 
with different species atomic ions[2ni El] which has been implemented for a range of 
masses 1/3 < ma/ < 2. Outside this range, ground-state-coohng becomes less 
effective, but might still be feasible. Suitable atomic ions range from 9 atomic mass 
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units (^Be"*") to 176 atomic mass units (^'^^Yb'^). Only one of the normal modes is 
required to be in the ground state for what follows. At this point the molecule is ideally 
in the electronic, vibrational and translational ground states, leaving only the rotational 
and, if the molecule has hyperfine structure, the hyperfine states undetermined. 

3. Rotational spectroscopy and cooling 

Manipulation of rotational states, ideally without altering the vibrational and the 
electronic state, might proceed by combining Raman-transitions p2j driven by optical 
frequency combs p3] with a technique similar to what has been used for state preparation 
and read-out in two- ion quantum logic clocks [21]. In this way, entropy can be removed 
from the molecular rotational degree of freedom by indirectly coupling it to the near 
perfect reservoir of a light field. This can be accomplished via closed-cycle electronic 
transitions of the co-trapped atomic ion. Raman transitions are a near-universal tool 
in molecular spectroscopy, because they can be driven with light sources that are not 
specific to the molecule under study. Direct electronic transitions in molecular ions 
typically require wavelengths shorter than 650 nm[2l], below the center wavelengths of 
common femtosecond laser frequency combs (mostly 800 nm or 1064 nm)[23]- Two comb 
beams can be produced either by splitting the output of one laser and frequency shifting 
one part[2S], or by using two separate combs with equal mode spacing and scanning the 
carrier-envelope phase of one or both. If the combs span about 30 THz, they are much 
broader than the transition frequencies between the highest lying thermally occupied 
rotational states that need to be addressed (typically below 300 K ks/h ^ 6.25 THz 
where ks and h are Boltzmann's and Planck's constant respectively[2l]). The two 
comb fields should counter-propagate along the ion crystal axis for maximal momentum 
transfer in transitions that change the state of one normal mode of the ion crystal 
while simultaneously driving a rotational state changing transition in the molecule (see 
Figjl]). Analogous processes have been utilized in experiments with a laser comb and two 
identical atomic ions where an approximately 12.7 GHz hyperfine Raman transition was 
driven simultaneously and coherently with the ion crystal mot ion [25]. Raman transitions 
of the molecule require that the resonance condition 

Num + = r^j,j±i,2 + nu± (4) 

is fulfilled. Here N and n are integers, is the mode spacing of both comb fields, 
set to be identical, and AcUo is the frequency shift between the two comb fields due to 
their individual carrier-envelope phases (or due to the frequency shift of one part of the 
output). The frequency difference between two rotational levels with quantum numbers 
J and J ± 1, 2 is denoted f2j,j±i,2 (for linear and spherical-top- molecules, J — j- J ± 1 is 
forbidden by Raman selection rules) and u± is either one of the normal mode frequencies 
of the two-ion crystal. Transitions with n = ("carrier") will have the largest transition 
rate, while the n-th "sideband" will be suppressed by the n-th power of the Lamb-Dicke 
parameter. For Raman-transitions induced by two counter-propagating beams with 
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Figure 1. Schematic setup of a quantum-limited molecular spectrometer. A known 
atomic ion (blue) and a, possibly unknown, molecular ion (red and blue stick model) 
are confined in the same harmonic trapping well. The atomic ion is initialized to a 
known internal state. Due to the Coulomb interaction between atomic and molecular 
ions the normal modes of vibration of the entire crystal can be sympathetically cooled 
to their ground state by coupling exclusively to the atomic ion with laser interactions 
(blue arrow). The rotational states of the molecule are then manipulated by inducing 
Raman-transitions due to the light fields of two counter-propagating frequency combs 
(orange arrows) of equal repetition rate aJm/(27r) and variable relative carrier-envelope- 
phase Awq. The photon recoils produced by Raman-transitions of the molecule can 
excite the shared crystal motion, and this excitation can subsequently be detected on 
the atomic ion (see text). Detection of crystal excitation thus establishes that the 
molecule has undergone a transition. 



wave- vector modulus fc/, the Lamb-Dicke parameters of the molecular ion are given 



At mode frequencies commonly used in resolved sideband cooling[T6], the ground state 
wavefunction extension of the molecule in either normal mode (typically a few nanome- 
ters) is much smaller than the light wavelength A, so higher-order sidebands are largely 
suppressed. The first sidebands (n = ±1) are the strongest motion-altering interaction 
and change the motional state by one quantum. 




(5) 
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For a very general rough estimate of achievable Raman rates, we can assume that 
the intermediate level in the Raman transition has a dipole matrix element of e ■ ao ~ 
8.48 10~^° C-m, where e is the electron charge and oq the Bohr radius. The transition 
rate on the n = carrier is then approximately [26] 

with A; the detuning of the comb center frequency to the resonance frequency from the 
molecular ground state to the intermediate level and i?av the average electric field per 
comb, set to be identical for both combs for simplicity. If each comb emits 1 W on 
average into a beamwaist of 20 /im around 800 nm and the dipole transition resonance 
is at 400 nm, the n = carrier Raman transition rate is |i7o|/(27r) ^ 824 kHz. For 
the n = ±1 sideband transitions this frequency must be multiplied by the Lamb-Dicke 
parameter 77, which typically ranges between 0.03 and 0.3, depending on the normal 
mode structure and the comb wavelength [see Eq.(|5])]. 



3.1. Rotational cooling for known molecular constants 

For known rotational constants and transition rates of the molecule, we can set the 
relative detuning of the combs such that two rotational levels with high thermal 
occupation probability are connected by a n = —1 resonance condition in Eq.Q. 
As long as the motion of the two-ion-crystal is in the ground state, this transition 
can proceed only one way, from the J to the J — 1, 2 states, while one quantum 
of motion is added to the mode at uj± (see Fig. [2]^a)), the reverse direction cannot 
conserve energy. Each attempt to excite a Raman transition between rotational states 
of the molecule with the comb-lasers would be followed by an attempt to drive a 
n — )■ (n — 1) resolved sideband on the cooling and detection ion (see Fig. [2|^b)). As 
long as the selected normal mode of crystal motion is in the ground state (n = 0), 
this transition is energetically forbidden, but if the preceding transition of the molecule 
has deposited a quantum of motion (n = 1), the electronic state of the atomic ion 
can be changed while this quantum is removed [21]. The electronic state-change of 
the atom can be detected by state-selective fluorescence on the cooling and detection 
transition, in analogy to quantum-logic ion clocks [Ml [21] (see Fig. gc)). Crucially, if 
a state change of the atomic ion is detected, this near-ideal measurement projects the 
rotational state of the molecule into a pure and known quantum state, reducing the 
entropy of the molecular state considerably in one successful attempt [27] . In case the 
excitation attempt is not successful, the addressed rotational state is almost certainly 
not occupied, so we can proceed to probing another highly occupied rotational state in an 
updated distribution that excludes the unsuccessfully probed level(s). Once the initial 
distribution of rotational states is projected into a pure state in this manner, we can 
drive the system straight into the rotational ground state with carrier-type transitions 
(n = in Eq.dl])), if the Raman transition rates are known. It might also be possible 
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Figure 2. Molecular rotational Raman transition and subsequent detection (energy 
differences not to scale), (a) The molecule undergoes a Raman-transition on the 
n = — 1 sideband of the selected normal mode of the two- ion crystal in Eq.Q. Since 
this mode is in the ground state, only the rotational transition J — > J— 1 is energetically 
allowed (solid lines) while the opposite direction is not resonant with an existing state 
(dashed lines), (b) The conditionally added quantum of motion is transferred into 
an excitation of the co-trapped atom by use of a (g,n = 1) — >■ (e,n = 0) resolved 
sideband transition, where g and e refer to two long-lived states in the atomic ion. 
(c) The electronic state of the atomic ion is detected by laser induced fluorescence 
on a cycling transition connecting g with the fast-decaying level /[Ml [H]. If no 
fluorescence is detected, the atom is found to be in excited state e. This near-ideal 
measurement projects the state of the molecule into J — 1. The atom can then be 
reset to g by optical pumping and the crystal can be re-cooled to the motional ground 
state. This resets atom and two-ion-crystal motion to the initial condition before step 
(a) without affecting the rovibrational or electronic state of the molecule. The reset 
removes entropy from the system and prepares it for further manipulations. 



to use techniques similar to stimulated rapid adiabatic passage (STIRAP)[T3l HU |28] 
to provide more robust carrier and sideband transitions without exact knowledge of the 
transition rates. 

3.2. Determination of transition rates 

If the Raman transition rates of the molecule are not known, they can be determined 
as follows: Once the molecular ion has been projected into rotational state J, we 
can probe the molecule on the {J,n = 0) — ?■ (J — l,2,n = 1) transition for several 
different excitation periods and transfer any motional excitation of the two-ion crystal 
into an electronic state change of the atom by inducing a red sideband transition of the 



atom, as in the preceding step (section 3.1) that narrowed the thermal J-distribution 



to just one state. If the atomic ion changes state in such an excitation attempt, the 
transition to {J — 1,2, n = 1) has most likely been made; therefore, the molecule is 
now in state J — 1,2 with high probability, and can be returned to J by re-cooling 
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the two-ion crystal to the n = ground state, followed by repeated attempts to do 
a ( J — 1, 2, n = 0) — 7- {J,n = 1) sideband transition (detuned by 2 ijj± to the blue 
from the transition previously attempted). By recording transition probabilities as a 
function of length of excitation in many repetitions of this cycle, the transition rate and 
subsequently the magnitude of transition matrix elements can be inferred. 



3. 3. Quantum-non-demolition measurement of the rotational state 

For known transition rates, one should be able to repeat the cycle of successive n = ±1 



transitions and detection of the molecular state-changes outlined in section 3^ several 
times. This ideally constitutes a quantum-non-demolition (QND) readout as we should 
leave a motional quantum that can be subsequently read out through the state of the 
atom in each cycle step. Every time this happens, we gain more certainty about the 
state the molecule is in, even if single measurements have only modest fidelity. An 
analogous state-detection scheme with atomic ions has reached a detection fidelity of 
99.94 %, with a fidelity of about 80 % for a single measurement [29]. Several cycles of 
such QND readouts also prepare the molecule rotation into a pure quantum state with 
higher confidence than that for a single detection of a transition. 



3.4- Determination of molecular constants 

If the rotational levels of the molecular ion are unknown, we can slowly scan the fre- 
quency difference Auo of the two combs over the comb spacing Um- At some point, a 
n = ±1 resonance will lead to a molecular Raman-transition that will leave a quantum 
in the selected mode of motion of the crystal. Detecting a state change in the atom 
after converting the motion into atomic excitation implies that the molecule is in either 
one of two J-levels that the previous Raman-interaction has resonantly connected, so 
the distribution of rotational states is narrowed to two levels. We can reset the atom 
state by optical pumping, and recool the crystal motion to the ground state, and then 
attempt to drive a molecular transition at a relative detuning that is 2 ll!± red-detuned 
from the previous pulse. Again, the transition is possible only if a quantum of motion 
is added that can be subsequently detected with the atom. Detection of the atomic 
state change would establish that the molecule has undergone a transition from the 
higher energy J to the J — 1 or J — 2 state (for linear and spherical-top molecules, 
the J — 2 state has been prepared with certainty). If several attempts to drive at this 
red-detuned frequency fail, the molecule is most likely already in the lower-energy state 
and the failed attempts were actually on a n = —3 sideband transition that is also 
energetically forbidden if the selected mode of crystal motion is in the ground state. In 
this case, attempts to drive transitions 2 ijj± blue detuned from the pulse that produced 
the original transition should excite the molecule to the higher-energy J state, while a 
motional quantum is deposited (n = 1) that can subsequently be converted to a state 
change of the atom and detected. In either case, the unknown frequency f2j,j±i,2 can 
be found from the mean of the two resonant frequencies in the preceding steps. It is 
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then possible to find the rate of these transitions and gain information about the related 
transition matrix elements as described in section 3.2 Further rotational levels may be 
identified by completing the scan of AcUo over the full comb mode spacing Um- After 
several distinct transitions are found, it should be possible to assign J and determine all 
rotational constants and centrifugal distortion constants. This can be done unambigu- 
ously, if the absolute frequency of every transition is known. This is possible even if N 
in Eq.Q is initially unknown. From the range of rotational transitions and assuming 
that Um is approximately 2tt x 1 GHz, can take values between 1 and 10^. After 
finding a resonance, we can change Um by a small amount of (say) order 10~^Um, easily 
observable by monitoring the repetition rate of the comb(s) on a fast photodiode. Then 
Auo can be scanned until the resonance condition is re-established. The closest integer 
to the fraction of the changes in Um and Auo is A^. We can eventually infer the full 
rotational spectrum of the molecule by systematically scanning the complete range Um 
and then determining A^ for each resonance found. With this knowledge, the spectrum, 
originally compressed into one range Um, can be deconvoluted to yield rotational states 
information with a resolution limited only by the lifetimes of the involved states. 



For simplicity, the preceding discussion has neglected any hyperfine structure that 
the molecule might possess. Hyperfine structure complicates the spectrum-finding 
procedure by adding more resonances and one more degree of freedom that needs to 
be determined. Still, it seems straightforward to incorporate hyperfine structure with 
slight modifications to the outlined procedures. 



4. MgH+ co-trapped with Mg+ 

Estimated parameters for the concrete case of a MgH"*" molecule, co-trapped with a Mg"*" 
atom establish the feasibility of the proposed scheme. The rotational state of MgH"*" has 
been previously prepared close to its ground state in an ensemble of molecular ions with 
dedicated laser systems and a destructive detection method [5]; therefore all relevant 
molecular parameters are reasonably well known. As a linear molecule, MgH"*" admits 
only J — 7- J ± 2 rotational Raman transitions. The low mass and small size produce 
a rotational constant of approximately 190 GHz and a centrifugal correction of order 

10 MHz. In thermal equilibrium with a room-temperature environment only the first 

11 rotational states have occupation probabilities above 1 %. The highest occupation 
probability is approximately 14 % for J = 4. The J — i- J±2 selection rule dictates that 
we should first try to prepare a state with even J to be able to connect to the ground 
state. The dipole transition matrix element between the ground X^S"*" and the lowest 
electronically excited A^S"*" state with an energy difference corresponding to 280 nm is 



about 1.57 e ao[30]. For the same laser parameters as in section 3.1, the carrier Raman 
transition rate, Eq.(|6]), is approximately 1.1 MHz. For a crystal consisting of MgH"*" and 
Mg"*" at a single Mg"*" trap frequency of =1 MHz the two normal- mode frequencies, 
Eq. g, are cu+ ^ 0.99 MHz and u_ ^ 1.72 MHz, with Lamb-Dicke factors, Eq. Q, of 
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77+ ~ 0.16 and //_ ~ 0.13 respectively. These numbers imply n = ±1 sideband transition 
rates of 174 kHz (complete rotation between states achieved in 1.43 /is) and 137 kHz 
(complete rotation in 1.82 /is), respectively. Together with initial ground-state cooling 
and subsequent transfer of motional excitation and state readout of the Mg"*" ion, one 
experimental cycle should take no more than 1 ms[TB]. An initial scan could involve 
100 attempts per comb detuning Acoo, yielding a scan speed of 10 frequency points per 
second. A 10 /xs window with about 10000 Raman pulses applied at 1 GHz repetition 
rate to the molecule would correspond to a Fourier limited width of about 10 kHz for 
each resonance; so to collect a reasonably densely sampled spectrum over 1 GHz, we 
would need 10^ points taken in 10^ seconds, or about 3 hours. This duration represents 
a worst-case estimate, where nothing about the rotational constants is known, and is 
at the same time compatible with achievable lifetimes in ion traps. If we assume that 
10 transitions with initial state population larger than 1 % are present, it should take 
less than 1 hour to identify the first resonance. With each newly identified resonance, 
finding the remaining resonances should take fewer attempts. In this fashion, it should 
be possible to determine a complete rotational spectrum in a few hours. Once the most 
populated sub-levels are determined, these could be directly prepared within a time 
span comparable to the rotational re-thermalization period of the molecule by repeat- 
ing excitation attempts until a transition is successfully made. An experiment in a low 
temperature environment would reduce the initial spread of J-levels considerably, while 
suppressing background gas collisions with the crystal and thus further increasing the 
ion lifetimes in the trap. 



5. Conclusions 

In the recent past, quantum control of neutral and ionized atoms has proven enormously 
fruitful. Controlling the quantum state of a multitude of molecules by use of the methods 
proposed here and also other methods that have recently emerged promises to take us to 
the next level of complexity. We can hope to keep and study the same molecule with a 
resolution limited only by natural linewidths under full control of its quantum states. It 
should also be possible to efficiently prepare certain rotational states, including the ro- 
tational ground state and coherent superpositions of these states. This may be further 
aided by a higher fidelity, quantum-non-demolition way of measurement /preparation 



that is described in section 3.3 After preparation, the molecular quantum state can be 
coherently manipulated with the same or other laser sources within the limits imposed 
by decay and re-thermalization of the involved states. For rotational states within 
the electronic and vibrational ground state manifold, re-thermalization is often on a 
timescale of many seconds. 



Aided by these tools, we might test quantum mechanics and its application to cou- 
pled systems to unprecedented levels and search for energy shifts due to parity violation. 
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a finite value of tfie electron electric dipole moment, or a variation of fundamental con- 
stants, such as the proton-to-electron mass ratio. It might also become feasible to induce 
highly controlled molecular collisions. Because reaction cross-sections can strongly de- 
pend on the vibrational and rotational quantum state, we might glean better control 
and understanding of chemical reactions in this way. 

During completion of this manuscript I became aware of similar work in the group 
of D. Matsukevich at the National University of Singapore [31]. 
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